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ABSTRACT 


Stability limits were determined, experimentally, for a 2 -inch (5. 08 cm) diameter 
rocket engine. Liquid oxygen and gaseous hydrogen were the propellants. Chamber 
pressures of 650 and 300 psia (4. 48xl0 6 and 2. 07xl0 6 N/m 2 ) were run with an oxidant- 
fuel ratio of 5. 0 and a characteristic length of 95 inches (2. 41). For comparison with 
experimental data, boundaries were generated on the analog computer using the chugging 
model proposed by the authors and calculated values of delay. Results agreed with 
regard to observed frequencies and boundary shape. Discrepancies in boundary location 
(required pressure drops) were attributed to a high combustion noise level. 


EXPERIMENTAL VERIFICATION OF A DOUBLE DEAD-TIME MODEL 
DESCRIBING CHUGGING IN LIQUID BIPROPELLANT 
ROCKET ENGINES 

by John R. Szuch and Leon M. Wenzel 
Lewis Research Center 

SUMMARY 

Low frequency combustion instabilities in liquid propellant rocket engines, commonly 
referred to as chugging, have been the subject of many analyses during the past two 
decades. The purpose of this investigation was to verify, experimentally, the chugging 
model previously proposed by the authors. The proposed double dead-time model con- 
siders each propellant to be acted upon by a discrete time delay between injection and 
burning in the combustion chamber. This differs from the commonly used single delay 
model where both propellants are acted upon by a single delay. 

Use of the double dead-time model results in a stability boundary having several dis- 
tinctive characteristics. More than one range of chugging frequencies may be observed 
along the boundary. Also, a reversal in slope of the boundary occurs at high fuel injector 
pressure drops, indicating the possibility of stabilization by decreasing the fuel injector 
pressure drop. 

Stability boundaries were determined, experimentally, for a 2-inch (5. 08 cm) diam- 
eter rocket engine. Liquid oxygen and gaseous hydrogen were the propellants. Bound 

aries were determined for two configurations: a chamber pressure of 650 psia 
fi 2 ^ 

(4. 48x10 N/m ), contraction ratio of 16. 8, and a chamber pressure of 300 psia 
(2. 07x10® N/m^), contraction ratio of 8. 4. The oxidant-fuel ratio was maintained 
around 5. 0 with a characteristic length of 95 inches (2. 41 m). 

For comparison with experimental data, boundaries were generated on the analog 
computer using the double dead-time model and calculated values of delay. Experimental 
and computer data agreed with regard to the observed frequencies (both higher and lower 
ranges were observed with each configuration) and the boundary shape. (The reversal in 
slope was observed. ) Discrepancies in boundary location (injector pressure drops 
required for stability) were attributed to a high combustion noise level. 



INTRODUCTION 


Low frequency combustion instabilities in liquid propellant rocket engines, commonly 
referred to as chugging, have been the subject of many analyses during the past two 
decades. The purpose of this investigation was to verify, experimentally, the chugging 
model advanced in reference 1. A comparison of the proposed double dead-time model 
and the commonly used single -delay model is given in figure 1. For the single -delay 
model, the injector flow rates are acted upon by a single delay, usually assumed to be 
made up of the governing vaporization time, mixing, and reaction times. For an oxidizer- 
vaporization limited system, this approach would be valid for the case of high fuel injec- 
tor pressure drop. Under these conditions, fuel flow rate would be relatively constant 
during chugging, approximating the monopropellant case. To analyze chugging over a 
wide range of operating conditions and especially when one of the propellants is introduced 
as a gas at low pressures, the single -delay treatment will not suffice. The double dead- 
time model considers each propellant to be acted upon by a discrete time delay between 
injection and burning in the combustion chamber. 

Figure 2 shows typical stability boundaries for both models, plotted as ratios of 
injector pressure drop to chamber pressure (AP/P). All symbols are defined in 
appendix A. Regions above and to the right of the boundaries correspond to stable oper- 
ation. For the case shown, the fuel is introduced as a gas with no vaporization time 
(r v ^ = 0). The single-delay boundary has the usual hyperbolic shape with the entire 
boundary corresponding to one chugging frequency. The double dead -time boundary has 
several distinctive characteristics. As the fuel AP/P is increased, the chugging fre- 
quency approaches a value (160 Hz) determined by the total oxidizer delay (2.75 msec) 
and the gas residence time. As the oxidizer AP/P is increased, the frequency 
approaches a value of about 360 hertz, determined by the total fuel delay (1. 0 msec) and 
the gas residence time. The boundary is made up of two portions, each corresponding 
to a different range of chugging frequencies. Multi -mode oscillations have been observed 
and reported in reference 2. The lower frequency portion of the boundary has a charac- 
teristic reversal in slope indicating the possibility of stabilizing an engine by decreasing 
the fuel injector AP/P over a limited range. Such behavior has been observed and 
reported in reference 3. 

The analysis of the bipropellant system considered in reference 1 was simplified by 
assuming constant injector cavity pressures (decoupled feed system). The system equa- 
tions were linearized by assuming small perturbations in the system variables about 
steady-state operating values. Efforts were made to achieve these conditions during the 
experiments reported herein. 

Experiments were conducted in the Rocket Combustion Laboratory at the NASA- 
Lewis Research Center. A 2-inch (5. 08-cm) diameter rocket engine was run at nominal 
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chamber pressures of 650 and 300 pounds per square inch (4. 48xl0 6 and 2. 07xl0 6 N/m 2 ) 
absolute. The propellants were gaseous hydrogen at room temperature and liquid oxygen 
(LOX) at 140° R (77. 8° K). The oxidizer -fuel mixture ratio was maintained around 5. 0 
with a total propellant mass flow rate of about 0. 7 pound per second (0. 317 kg/sec). 
Transitions from stable to unstable operation were accomplished by varying the injector 
geometries while maintaining the primary system variables (chamber pressure, flow 
rates) constant. Stability boundaries were determined for both configurations and com- 
pared with boundaries established on the analog computer using the double dead-time 
model. 


APPARATUS 

The flow system is shown schematically in figure 3. Pressure-regulated propel- 
lants, at ambient temperatures, were supplied to the system from storage bottles. The 
gaseous oxygen was liquefied in a liquid nitrogen bath. Valve and ignition sequencing 
was controlled by a preset electric timer. The chronology of ignition was as follows: 

(1) spark plug ignition and hydrogen bleed flow, (2) main oxidizer flow, and (3) main 
hydrogen flow. Relative timing between events was adjusted to give a smooth start. 


Injector 

A cut-away drawing of the injector is shown in figure 4. To conform with the model 
assumed in reference 1, it was necessary to maintain constant injector cavity pressures. 
That is, the pressure perturbations induced in the injector cavities during chugging had 
to be sufficiently attenuated. This attenuation was accomplished on the hydrogen side by 
adding a large volume (300 in. 4. 92xl0 -2 m 2 ) immediately upstream of the hydrogen 
injector manifold (see fig. 3). This method was not practical for the oxygen side, how- 
ever, due to the high bulk modulus of liquid oxygen. Attenuation of the oxygen injector 
perturbations was provided by a thin (0. 063 in. ; 0. 16 cm) inconel diaphragm (A). It was 
not possible to design a diaphragm with sufficient flexibility and yet sufficient strength to 
contain full manifold pressure. Therefore, most of the manifold pressure was balanced 
by gas pressure in volume (B). Volume (B) was pressurized, prior to running, to the 
anticipated manifold pressure minus about 15 psi (1. 04xl0 5 N/m 2 ). The thick perforated 
plate (C) supported the diaphragm against this preload. Careful design was required to 
insure that the diaphragm compliance would not couple with the inertance of the holes in 
plate (C) to produce a resonance in the frequency range of interest (30 to 300 Hz). 
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The hydrogen injector pressure drop was manipulated by changing the thickness 
(hence, porosity) of the transpiration-cooled faceplate (D). Thicknesses ranging from 
0. 125 to 0. 250 inch (0. 318 to 0. 635 cm) were used to give the desired range of injector 
pressure drops. 

The orifice plate (E) was used to vary the oxygen injector pressure drop. Flow 
diameters from 0. 076 to 0. 1285 inch (0. 193 to 0. 326 cm) were used. Minimum pressure 
drop was achieved with no orifice in the system. The injection tube (F) was made about 
ten diameters long in an attempt to maintain constant liquid jet properties at the chamber 
entrance. 


Engine 

An uncooled combustion chamber was used. Chamber wall thickness was 1. 0 inch 
(2. 54 cm). The inside diameter was 2. 0 inches (5. 08 cm) and the length was variable 
from 2. 1 to 12. 5 inches (5. 33 to 31. 8 cm). By changing chamber length, gas residence 
time (for stability studies) and burning length (for performance studies) could be eval- 
uated. 

With constant flow rates, chamber pressure was varied by changing the throat diam- 
eter of the water-cooled exhaust nozzle. Data were taken with throat diameters of 
0. 487 and 0.689 inch (1.24 and 1.75 cm), yielding nominal chamber pressures of 650 
and 300 psia (4. 48xl0 6 and 2. 07xl0 6 N/m 2 ), respectively. 

Instrumentation 

The LOX flow rate was measured by means of a turbine-type flowmeter located in the 
nitrogen bath (see fig. 3). The gaseous hydrogen flow rate was measured with a critical 
flow venturi, located between the hydrogen fire valve and the 300 cubic inches 

O «} 

(4. 92x10 m ) accumulator. Steady-state pressures were measured with strain gauge 
type transducers. Dynamic pressure measurements were made in the chamber and both 
injector cavities using piezoelectric transducers. The frequency responses of all 
dynamic pressure transducers were flat well beyond the frequency limit of this study 
(300 Hz). Output signals from all steady -state instrumentation were displayed on 
electrically -clamped meters. All transducer outputs were recorded on oscillograph 
paper. In addition, the output signal from the dynamic chamber pressure transducer was 
recorded on an FM tape recorder. The recorded signal was played back through a fre- 
quency analyzer, consisting of 30 filters with bandwidths of 1/3 octave and center fre- 
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quencies ranging from 40 to 31 500 hertz. Maximum steepness of the filter frequency 
characteristic was 120 decibels per octave. 


PROCEDURE 
Cold Flow Tests 

Cold flow tests were conducted prior to the taking of stability data. With no hydrogen 
flow and the LOX diaphragm supported by high gas pressure, oxygen was flowed through 
the injector. LOX injector pressure drop against LOX flow rate is plotted in figure 5 for 
orifice diameters between 0. 076 and 0. 1275 inch (0. 193 and 0. 326 cm). As expected, 
pressure drop was proportional to the square of flowrate for a given orifice diameter. 

With no LOX flow, room temperature hydrogen was exhausted through the fuel injec- 
tor to the atmosphere. Faceplate thicknesses of 0. 141 and 0. 175 inch (0. 358 and 
0. 444 cm) were run. Due to the compressibility of the hydrogen, flow through the face- 
plate is a function both pressure drop and downstream conditions. Assuming that 
(1) hydrogen temperature remains constant, and (2) pressure is distributed linearly 
through the faceplate, hydrogen flow rate can be calculated as follows: 


w. 


= < C d A ) R ( 2 Vf,av AP tt ) 1/2 

/ 2g X 1 / 2 / v 

= (C ,A) — \ (P AP Tf ] 

d R \R f T fy ) \ f > av 

= K R [< P K - P c>< P H + P c>] 1/2 


1/2 


( 1 ) 


where 

p p 

g c force-mass conversion factor, 396 (lb m )(in. )/ (lb f )(sec ); 1 (kg)(m)/(N)(sec ) 

Wj hydrogen mass flow rate through the Rigimesh faceplate, lb m /sec; kg/sec 

K r faceplate coefficient (function of thickness), (lb m )(in. 2 )/ (sec)(lb^); (kg)(m 2 )/(sec)(N) 

9 p 

P H fuel injector pressure, lb^/in. ; N/m 

2 2 

P c chamber pressure, lb^/in. ; N/m 
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Figure 6 contains a plot of hydrogen flow rate against the quantity [(P Tf - P ) 
i j/ 2 ^ ^ c 

(Pjj + P c )J ' for room temperature hydrogen exhausted to the atmosphere through 

specified faceplate thicknesses. The linearity of these plots and equation (1) permit 

computation of expected injector pressure drops for selected faceplate thicknesses and 

chamber pressures of 650 and 300 psia (4. 48x10® and 2. 07x10® N/m^). 

Figure 7 contains resultant information for room temperature hydrogen flowing at 

0. 11 pound per second (0. 0499 kg/ sec). Also shown are actual data points obtained 

during the hot run phase of the program. 


System Definition 

In an attempt to reduce the number of variables governing the vaporization process, 

it was decided to maintain constant LOX flow rate (hence, injection velocity), mixture 

ratio, and characteristic length, L* while running at each of the two throat sizes. With 

a throat diameter of 0. 487 inch (1. 24 cm), a chamber length of 4. 1 inches (10. 4 cm) 

was selected. This configuration yields an L* of 91. 0 inches (2. 31 m) and will be 

referred to as configuration A. With this configuration, a LOX flow rate of 0. 55 pound 

per second (0.249 kg/sec) was required to reach a steady-state chamber pressure of 

650 psia (4. 48x10° N/m ) at a nominal mixture ratio of 5. 0. 

With the available chamber lengths and a throat diameter of 0. 689 inch (1. 75 cm), 

the closest L* that could be run was 95. 5 inches (2. 43 m) with a chamber length of 

9. 85 inches (25. 0 m). Running this configuration at a LOX flow rate of 0. 55 pound per 

second (0. 249 kg/ sec) and a mixture ratio of 5. 0 resulted in a nominal chamber pressure 
6 2 

of 300 psia (2. 07x10° N/m;. This configuration will be referred to as configuation B. 
Further description of both configurations may be found in table I. 


Stability Investigation 

For both configurations, orifice diameters and faceplate thicknesses were selected 

using figures 5 and 7, respectively, to give a range of injector pressure drops from 

0.05 P to 1. 0 P . In the case of the LOX injector however a minimum of about 
c 5 c 2 

50 psi (3. 45x10° N/mj could be obtained at the prescribed flow rate by removing the 
orifice. 

Since system operation was open-loop, it was necessary to iterate to the desired 
steady-state operating conditions. For a given system configuration, initial tank pres- 
sures were selected and set. With the ignition sequence controlled by the electric timer, 
the engine was fired for approximately 4 seconds. Timing was adjusted to give 
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slightly more than 2 seconds of steady -state operation. Adjustments were made in the 
tank pressures so as to decrease the errors between measured and desired flow rates. 

In general, LOX and fuel mass flow rates were held within 10 percent of 0. 55 and 0. 11 
pound per second (0. 249 and 0. 0499 kg/sec), respectively. Destabilization of the engine 
was accomplished by increasing the LOX orifice size while maintaining faceplate thick- 
ness and the other system variables constant. 

The electric timer was also used to gate the output of the dynamic chamber pressure 
transducer to the tape recorder. The gating circuit also provided 2 seconds of refer- 
ence ground both before and after the steady-state data. Figure 8 shows a typical start 
transient, recorded on oscillograph paper with the timer events noted. 


Performance Investigation 

A model describing the vaporization process for a single droplet has been developed 
by Priem and Heidmann (ref. 4). For the case of a completely vaporized fuel, combus- 
tion efficiency is completely defined by the percentage of the oxidizer mass vaporized 
within the chamber. Heidmann and Wieber (ref. 5) have been able to correlate frequency 
response data obtained with this model with that obtained using a pure delay model with 
delay equal to the time required to vaporize 50 percent of the mass of an injected droplet. 
From these results, it follows that (1) at a chamber pressure of 300 psia (2. 07x10 
N/m^), (2) a mixture ratio of 5. 0, and (3) with a chamber length yielding a combustion 
efficiency of 62 percent, 50 percent of the oxidizer mass will be vaporized within the 
chamber. To determine this length, of the oxidizer mass will be vaporized within the 
chamber. To determine this length, a faceplate thickness of 0. 165 inch (0. 419 cm) and 
an orifice diameter of 0. 098 inch (0. 249 cm) were selected to give stable operation with 
a throat diameter of 0. 689 inch (1. 75 cm). Chamber lengths, ranging from 2. 1 to 12. 5 
inches (5. 33 to 31. 8 cm) were run with tank pressures adjusted to give a chamber pres- 
sure of 300 psia (2. 07x10® N/m^) and a mixture ratio of 5. 0. 


Data Analysis 

Frequency analysis of the recorded chamber pressure data showed that, in all cases, 
the chamber pressure oscillation could be considered as having two frequency compo- 
nents. Using the frequency analyzer filters, the two filtered signals for each run were 
rerecorded on oscillograph paper for further analysis. Figure 9 shows the envelope of 
a typical filtered signal observed during chugging. Due to the observed randomness in 
chamber pressure amplitude, a stability criterion based on a time -average over a set 
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portion of the run was selected. Any oscillation having an rms amplitude exceeding 
10 percent of the mean chamber pressure for 1 second was considered unstable. 

A total of 76 experimental runs was made with configuration A and 53 with configur- 
ation B. It was felt that some realistic criteria should be applied to the data to separate 
’’good” from ’’bad” data points. The following criteria were selected: (1) orifice 
repeatability, (2) faceplate repeatability, and (3) performance repeatability. Since cold - 
flow tests established the oxygen flow through the orifice as having a square law charac- 
teristic and since a constant density fluid was assumed, the following parameter was 
selected as a measure of orifice repeatability: 

C do <Po> 1/2 =W 2 So< P Io- P c>' 1/2 <2) 

Faceplate repeatability (for a given thickness) was based on the parameter as 

defined in equation (1). From observing the system behavior over a wide range of oper- 
ating conditions, combustion efficiency was found to be a function of the stability of the 
system. For a particular stability classification (lower or higher frequency, stable or 

unstable), combustion efficiency was selected as a measure of system repeatability. The 

1/2 -1 

parameters (p Q ) ' , K R , and ri c * were assumed to be normally distributed about 

some mean value with an allowable spread of ±1. 5 times the computed standard deviation. 
Based on the analysis, a total of 18 runs with configuration A and 10 runs with configur- 
ation B were judged to be "bad”. Pertinent data for the acceptable runs are listed in 
tables IE and HI. A summary of the data analysis is given in appendix B. System gains 
(injector AP/P) for both configurations were normalized to O/F = 5. 0 conditions as 
described in appendix C. 


RESULTS AND DISCUSSION 
Chugging Results 

Figure 10(a) shows the 58 acceptable runs for configuration A plotted on a AP/P 
map. A boundary, based on the 10 percent rms criterion, has been drawn between 
those points classified as stable or unstable. Corresponding chugging frequencies, for 
points near the boundary, have been noted on the figure. Oscillations in both the 70 and 
170 hertz ranges were observed. Note the characteristic reversal in slope of the lower 
frequency portion of the boundary. At an oxidizer AP/P of 0. 4, a reduction in the fuel 
AP/P from 0. 4 to 0. 2 resulted in a stabilization of the engine. 
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Figure 10(b) shows the 43 acceptable runs for configuration B. A boundary has been 
drawn based on the 10 percent rms criterion. Unlike the high pressure case, however, 
the boundary has an irregular shape. A bending back of the boundary takes place at both 
extremes of fuel injector pressure drop. Since control of fuel velocity was sacrificed 
for ease of changing the injector pressure drop, fuel velocity is not constant over the 
range of fuel injector AP/P’s. The resulting effect on oxidizer drop size and vapor- 
ization rate may be the cause of the irregular boundary shape. Another possible cause 
might be coupling of the feed system in spite of the efforts made to prevent such an 
occurrence. Chugging frequencies in both the 40 and 100 hertz ranges were observed, 
together with the reversal in slope of the low frequency portion of the boundary. 


Performance Results 

Figure 11 contains a plot of combustion efficiency against chamber length for a 

6 2 

chamber pressure of 300 psia (2. 07x10 N/m ). Results obtained in reference 6 sub- 
stantiate the low performance obtained with this injector. A combustion efficiency of 
75 percent would be expected with a chamber pressure of 300 psia (2. 07x10 N/m ), 
chamber diameter of 2. 0 inches (5. 08 cm), chamber length of 8. 0 inches (20. 3 cm), 
throat diameter of 0. 75 inch (1. 90 cm), mixture ratio of 5. 0, and a single -element, 
concentric tube injector. A slightly lower efficiency (71. 1 percent) was obtained with 
configuration B which approximates these conditions with the single -element, 
transpiration-cooled injector. 

Based on the vaporization model of Priem and Heidmann, the length required to 
vaporize 50 percent of the oxidizer droplet mass (« * = 0. 616) is 5. 5 inches (14. 0 cm). 

If one assumes that the average droplet velocity over this length is equal to the injection 
velocity (ref. 4), the vaporization time can be computed from l ^vr* = 48 milliseconds 

for configuration B. This result, when applied to the vaporization model of Priem and 
Heidmann, yields a vaporization time of 4. 62 milliseconds for configuration A. 

Values of mixing and reaction times were calculated for both configurations using the 
calculated values of vaporization time and the total oxidizer delays required to satisfy the 
phase requirement for neutral stability at high values of fuel AP/P and the observed 
chugging frequencies. Mixing times of 3. 8 and 2. 1 millisecond were determined for 
configurations A and B, respectively. 


Computer Simulation 

For comparison with experimental data, stability boundaries were generated using 
an analog computer simulation of the double dead-time model. The Analog Computer 
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Facility of the Lewis Research Center was used. Boundaries were obtained at chamber 
pressures of 650 and 300 psia (4. 48x10® and 2. 07x10® N/m^) with a LOX flow rate of 
0. 55 pound per second (0. 249 kg/sec) and a mixture ratio of 5. 0. The nonlinear differ- 
ential equations presented in reference 1 were programmed on the computer, using a 
1:1 time scale. A variable diode -function generator was used to represent the C^-O/F 
relation. Two variable delay lines were used to represent the total oxidizer and 
fuel delays. A one lump representation of the oxygen feed system (line, fire valve) was 
developed from measured line dimensions and valve calibrations. A summary of the 
equations, represented on the computer, is given in appendix D. 

A stability boundary was obtained on the computer by (1) selecting values for vapor- 
ization, mixing and gas residence time, (2) setting the value of K R to give the desired 
fuel injector pressure drop, (3) continuously varying K Q (LOX injector pressure drop) 
while adjusting tank pressure to maintain flow rate constant, (4) combining the chamber 
pressure signal with the output of a sinusoidal oscillator to produce a lissajous figure on 
an oscilloscope screen, (5) adjusting the oscillator frequency until a stable pattern 
appeared on the screen (the corresponding frequency being equal to the chugging fre- 
quency), (6) repeating the process at selected values of K R to form a boundary, 

(7) readjusting the delays and/or gas residence time to match the experimentally observed 
chugging frequencies over the entire boundary. 

Table IV gives the required values of delay together with the initial calculated values. 
LOX vaporization times were within 7 percent of those calculated using the vaporization 
model of Priem and Heidmann. All delay values were within 10 percent of the initial 
calculated values. 

In an attempt to match the randomness of the chamber pressure oscillations, the 
physical system was simulated with white noise superimposed on the products of combus- 
tion. The presence of a high combustion noise level could be attributed to the coarse 
injection pattern. The noise generator output was adjusted to give approximately 7 per- 
cent stable oscillations in chamber pressure with high injector pressure drops. This 
approximated the observed experimental noise level. Figure 12 gives a comparison 
between experimental and noise-induced chamber pressure oscillations at a fixed fuel 
injector pressure drop with variable LOX injector pressure drop. The 0. 25 level of 
oxidizer AP/P corresponds to the inherent stability boundary (no disturbance). It is 
seen that the transition from stable to unstable operation is a gradual one with amplitudes 
of oscillation ranging from 7 to 36 percent of the mean chamber pressure. Using the 
10 percent rms criterion for stability, a shift occurs in the apparent boundary position 
from 0. 25 to 0. 40. 

Figure 13 shows the sensitivity of the entire computer boundary to the input noise 
level for configuration A. It is seen that the introduction of noise into the system results 
in a shift in the entire boundary for a fixed stability criterion. A slight increase in the 
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boundary frequencies was noted with increased noise levels. 

Figure 14(a) gives a comparison between the experimental boundary for configur- 
ation A and the computer boundary obtained with the 7 percent noise level. Chugging 
frequencies between 63 and 85 hertz and between 147 and 210 hertz were observed both 
experimentally and on the computer. In reference 3, a ratio of 2. 7 between higher and 
lower frequencies has been noted during multi-mode oscillations. Data for that study 
were taken at high values of fuel AP/P so as to approximate a monopropellant system. 
An average ratio of 2. 5 was observed experimentally with configuration A at values of 
APjj/P c above 0. 5. Excellent agreement between the experimental and computer 
boundary shapes was obtained for configuration A. As predicted by the double dead-time 
model, a reversal in slope occurs for the lower frequency portion of the boundary with 
an intersection with the higher frequency portion of the boundary occurring at a fuel 
AP/P around 0. 15. Deviations between the boundaries at low fuel AP/P's may be due 
to the effect of fuel velocity on the oxidizer drop size and vaporization rate. The pos- 
sibility also exists that backflow through the fuel injector occurs at low values of fuel 
AP/P during chugging with a stabilizing effect. Excellent agreement inboundary posi- 
tion was obtained with the 7 percent noise level. 

Figure 14(b) gives a comparison between the experimental boundary for configur- 
ation B and the computer boundary obtained with the 7 percent noise level. Chugging 
frequencies between 39 and 46 hertz and between 99 and 112 hertz were observed. An 
average frequency ratio of 2. 5 was observed at values of fuel AP/P above 0. 5. The 
bending back of the higher frequency portion of the experimental boundary may be 
attributed to velocity effects or possible coupling of the feed system. A higher percent 
noise level would be required to match the boundary position indicating that the absolute 
noise level is a function of the injector configuration and not chamber pressure. A 
15 percent noise level with configuration B would be equivalent to the 7 percent noise 
level observed with configuration A. Both characteristics of the double dead-time model 
were exhibited by both boundaries however. A reversal in slope of the lower frequency 
portion of the boundary together with two ranges of frequency along the boundary were 
observed. 


SUMMARY OF RESULTS 

Stability limits were determined, experimentally, for a 2-inch (5. 08-cm) diameter 
rocket engine. Liquid oxygen and gaseous hydrogen were the propellants. Boundaries 
were determined for two engine configurations: a chamber pressure of 650 psia 
(4. 48x10 N/m ), contraction ratio of 16. 8, and a chamber pressure of 300 psia 

C Q 

(2. 07x10 N/m 4 ), contraction ratio of 8. 4. The oxidizer-fuel mixture ratio was main- 
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tained at 5. 0 with a characteristic length of about 95 inches (2. 41 m). For comparison 
with experimental data, stability boundaries were generated on the analog computer using 
the chugging model proposed by the authors in reference 1 and values of combustion delay 
determined from an existing vaporization model. The results of this study can be sum- 
marized as follows: 

1. The validity of the double dead-time model has been demonstrated. Both configur- 
ations investigated exhibited the characteristic behavior of the model. That is, more 
than one range of frequency along the stability boundary was observed together with a 
reversal in slope of the lower frequency portion of the boundary. 

2. Values of LOX vaporization times required to match the observed chugging fre- 
quencies were within 7 percent of those predicted by Priem and Heidmann. 

3. Mixing and reaction times must be inferred from the calculated values of vapor- 
ization time and the observed chugging frequencies. 

4. The choice of a single element injector configuration resulted in a high combus- 
tion noise level requiring the selection of an arbitrary time-average criterion for 
stability. Boundary positions were matched using the model and an experimentally 
determined 7 percent noise level. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 23, 1968, 

128-31-01-13-22. 
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APPENDIX A 


A cross-sectional area, in 2 ; m 2 

orifice coefficient 

C & theoretical characteristic veloc- 

ity, in/sec; m/sec 

d cross-sectional diameter, in; m 

f experimentally -observed fre- 

quency, Hz 

f' computer -observed frequency, 

Hz 

G system gain 

g c force-mass-conversion factor, 

386 (lb m )(in. )/(lbf)(sec 2 ); 

1 (kg)(m)/(N)(sec2) 

K Q oxidizer orifice coefficient, 

( lb m )( in )/( sec )( lb f) 1/2 ; 

0^)(m)/(sec)(N) 1 /2 

K r fuel faceplate coefficient, 

( lb m )( in - 2 )/( lb f)( sec ); 

(kg)(m 2 )/ (N)(sec) 

l length, in. ; m 

L* characteristic length, in. ; cm 

N number of data points in specific 

sample 

n white noise signal, lb m /sec; 

kg/sec 

O/F mixture ratio 

p 

P static pressure, lb f /in 

P time -derivative of pressure, 

lb f /(in. 2 )(sec); N/(m 2 )(sec) 


gas constant, in./°R; m/°K 

mixture ratio notation O/F 

root-mean-square amplitude of 
chamber pressure oscillation, 
lb f /in. 2 ; N/m 2 

estimated standard deviation of 
an assumed random parameter 
x, units of x 

s Laplace operator, sec"* 

T temperature, °R; °K 

t time, sec 

t R fuel faceplate thickness, in; m 

3 3 

V volume, in. ; m 

v velocity, in/sec; m/sec 

w mass flow rate, lb/sec; 

kg / sec 

Wj fuel mass flow rate, K^^P R - P 2 , 

lb m / sec; kg / sec 

w Q oxidizer mass flow rate, 

K o^Io"^ lb m /sec; kg/sec 
w time derivative of mass flow 

p p 

rate, lb m /sec ; kg/ sec 

x general system parameter x, 

units of x 

x mean value of parameter x, 

units of x 

a variable mixture ratio gain 

factor 


SYMBOLS 

R 

(ft 

rms 
S 
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AP 

pressure drop (assumed positive), 

i 

i measurements of variable 


lb^/in. N/m^ 


(i = 1, 2, . . . N) 

r 

specific heat ratio of combustion 

If 

fuel injector 


products 

Io 

oxidizer injector 

%* 

combustion efficiency 

3 

oxidizer jet 

p 

Q 

fluid mass density, lb m /in. ; 

LO 

oxidizer line 


kg/m 3 

m 

mixing 

T 

time delay, sec 

n 

exhaust nozzle 

6 s 

theoretical gas residence time, 

o 

oxidizer 

sec 

R 

transpiration-cooled faceplate 


Subscripts: 

t 

exhaust nozzle throat 



TO 

oxidizer tank 

av 

average 

vf 

fuel vaporization 

b 

products of combustion 

VO 

oxidizer vaporization 

c 

combustion chamber 

50 

50-percent mass vaporized 

f 

fuel 
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APPENDIX B 


DATA ANALYSIS 


x = selected parameter 
th 

Xj = i measurement of parameter (i = 1, 2, 3, . . . N) 
N = number of measurements 


x = mean value of parameter = 


N 

L x. 

i=l l 

N 


S x = estimated value of standard deviation of parameter 


N 

Z (x.-x) 2 
i=l 


-I 1/2 


(N - 1) 


Configuration A Analysis 


1. Orifice repeatability 


- w 172 = v; 1 K< p io- p c>]' 1/2 
1/2 


N = 76 


lb 


x = 0. 164 


m 


in. 


lb 


3/2’ 


27.3 


kg 


1/2 


3/2 


m 


1/2 


S = 0. 0212 -g 1 — ; 3.53 

Y « /«. 7 


kg 


1/2 


in 


3/2 


m 


3/2 


number of excluded points = 7. 
2. Faceplate repeatability 


x = w,- 1 (p2 . p2) 1/2 = K^ 1 
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Thickness 

N 

10‘ 5 

— 

X 

10~ 2 

S x 

Number of 
excluded 
points 

in. 

cm 

(lb^)(sec) 

(N)(sec) 

(lb f )(sec) 

(N)(sec) 




(lb m )(in. 2 ) 

(kg)(m 2 ) 

(lb m )(in. 2 ) 

(kg)(m 2 ) 


0. 125 

0.318 

6 

97.6 

14. 8 

7. 80 

1. 19 

0 

. 137 

.348 

5 

88.8 

13. 5 

3.89 

.591 

0 

. 141 

.358 

9 

67.8 

10.3 

2.99 

.454 

0 

. 148 

.376 

9 

47. 1 

7. 16 

1.92 

.292 

2 

. 154 

.391 

9 

37.2 

5.65 

1.42 

.216 

0 

. 158 

.401 

9 

34. 1 

5. 18 

2. 17 

.330 

0 

. 165 

.419 

9 

25. 4 

3. 86 

1. 04 

. 158 

1 

. 175 

.444 

7 

21. 0 

3. 19 

1.81 

.275 

0 

. 185 

. 470 

6 

17. 0 

2. 58 

1.82 

.277 

0 


3. Performance repeatability 

x = rj c * = A t g c P c (w o + Wf )" 1 C - 1 


Lower 

frequency 

classification 

Higher 

frequency 

classification 

N 

X 

S 

X 

Number of 
excluded 
points 

Stable 

Stable 

45 

0.739 

0. 0515 

6 

Unstable 

Unstable 

1 

.675 


0 

Stable 

Unstable 

4 

.628 

. 0660 

0 

Unstable 

Stable 

16 

.691 

. 0626 

2 


Configuration A total excluded points = 18. 


Configuration B Analysis 

1. Orifice repeatability 

* = C d<f 0 ) 1/2 = W 2 S C [< p Io - P c »]' 1/2 

N = 53 
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x = 0. 154 (ib^^n. _3/2 ); 25. 6 (kg 1/2 ) (m' 3/2 ) 

S x = 0. 0214 (lb^ 2 ) (in. _3//2 ); 3. 56 (kg l/2 )(m“ 3/2 ) 

number of excluded points = 5. 

2. Faceplate repeatability 


Thickness 

N 

icr 

in. 

cm 


(lb f )(sec) 




(lb m )(in. 2 ) 

0. 141 

0.358 

6 

63.0 

. 148 

.376 

7 

42.3 

. 154 

.391 

5 

35.3 

. 158 

.401 

5 

29.7 

. 165 

.419 

6 

22.8 

. 175 

.444 

6 

18.5 

. 185 

.470 

4 

14.9 

.200 

. 508 

2 

11.0 

.217 

.551 

3 

9.67 

.250 

.635 

4 

7.80 


3. Performance repeatability 

x = rj c * = A t i 


X 

icr 2 

S x 

Number of 

(N)(sec) 

(lb f )(sec) 

(N)(sec) 

excluded 

points 

(kg)(m 2 ) 

(lb m )(in. 2 ) 

(kg)(m 2 ) 


9. 58 

0. 693 

0. 105 

0 

6.43 

1.87 

.284 

1 

5.37 

.357 

.0543 

0 

4.51 

.430 

.0654 

0 

3.47 

.380 

.0578 

0 

2.81 

. 590 

. 0897 

0 

2.26 

1.05 

. 160 

0 

1.67 

.424 

.0644 

0 

1.47 

. 141 

. 0214 

0 

1. 19 

.687 

. 104 

t 

0 


c P c (”o + "'f)’ 1 c th 1 


Lower 

Higher 

N 

_ 

X 

S x 

Number of 

frequency 

frequency 




excluded 

classification 

classification 




points 

Stable 

Stable 

20 

0.711 

0. 0391 

2 

Unstable 

Unstable 

10 

.689 

.0286 

1 

Stable 

Unstable 

11 

.690 

.0302 

0 

Unstable 

Stable 

6 

.693 

.0347 

1 


Configuration B total excluded points = 10. 
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APPENDIX C 


NORMALIZATION OF ENGINE GAINS FOR OFF-MIXTURE RATIO CONDITIONS 


From the linearization of the system equations as described in reference 1, the 
characteristic equation for the closed-loop system can be written as: 


where 


0=1 + 


e" TmS L „ _T vo s , ^ 

I Gr e + G f 


(0g S + 1) 


G_ = 


0P dw 
c o 


° 0w„ d AP t (w + w,) 

c Io ' o 1 


1 + 


(tfl + 1) 9C th 


dw. 


C* 

C th 


3(fi. d AP 


Io 


G f = 


9P dw, 
c f 


0C. 


(R((R + 1) th 


dw. 


0(R d AP. 


0w f d AP h (w q + w f ) 
if w Q is proportional to ^Pj^^ 2 ; w f is proportional to (p^ - P 2 ^ ^ ; Pj Q , P 


nontime varying 


d »o _ *o 
d AP Io " 2 A Plo 

d w, w,P 

f_ f c 

d AP k AP H (AP If + 2P c ) 


Therefore, 


a 

G = — 

° AP. 


Io 


where 
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where 


_ (ft t (ft 9C th 


((ft + 1) 

C* 

°th 

3<ft 


Of 


APjc 
2 U 

P c 

/ ap h 

\ 2P c 

+ 1 ) 

1 

(ft 

9C ?h 

(<R + 1) 

Wh 

3<ft 


For the same system gains at off -mixture ratio conditions, the AP/P's can be 
normalized to (R = 5. 0 conditions by: 



AP n 

(i + iP “\ 

a f AP 

1 <R=5 If 

fi ♦ A M 

p 

\ 2P / 

«f p c 

1 (ft c 

V 2P / 

C 

' c Aft=5 

\ c / c 
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APPENDIX D 


SUMMARY OF EQUATIONS REPRESENTED ON THE ANALOG COMPUTER 


p 

1. Pressures in lb f /in. ; flow rates in lb /sec 


m 


P TO ” P LO ~ 155 w TO + 5# 90 ™TO 

P LO = 4 - 01Xl ° 6 f ^ TO - W LO> dt 

P LO ' P Io = 400 ^LO + 9 - 36 ^LO 
P lQ = 1. 37X10 3 f (w LQ - w Io ) dt 


P T - P = K w? + 0. 214 w Tr% 
Io c o Io Io 


e g p c + P c = xf- ho<‘ - T vo - T m> + "'rf< t - T m> + ”<*>] 
A t g c 

• / 2 2\ 4/ ^ 
w If = K R ( P H ” P c) 

Q/ F - *I° (t - T v° • T m> 

®H (t • T m> 


p 

2. Pressures in N/m ; flow rates in kg/sec 

p To - p LO = 5 - 19xl ° 6 + 8 - 97xl ° 4 w TQ 

P LO = 6> 10xl ° 10 / ^TO “ dt 

P LO - P Io = h 67Xl ° 7 ™LO + !• 42X1 ° 5 ^LO 
P Io = 2. 08x1 o' 7 /(<>, LQ - w lQ ) dt 

p Io - P c = Vlo + 3 - 25x11,3 *Io 
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TABLE I. - DESCRIPTION OF ENGINE CONFIGURATIONS STUDIED 


! 

Configuration 


A 

B 

Chamber diameter, d , in.; cm 

2. 00; 5. 08 

2. 00; 5. 08 

Throat diameter, d t , in. ; cm 

0. 487; 1.24 

0.689; 1.75 

Chamber length, l , in.; cm 

4. 10; 10.4 

9. 85; 25.0 

Characteristic length, L*, in. ; m 

91.0; 2.31 

95.5; 2.43 

Theoretical gas residence time, 



0g, at O/F = 5, msec 

2.40 

2.50 

Nominal mixture ratio, O/F 

5. 00 

5. 00 

Nominal chamber pressure, P , 



lb|/ in. N/ m^ 

650; 4. 48xl0 6 

300; 2. 07xl0 6 

Nominal LOX mass -flow rate, 



w o , lb m /sec; kg/sec 

0. 55; 0.249 

0.55; 0.249 


CO 


TABLE H. - CONFIGURATION A DATA 


Run 

Chamber pressure, 
P c 

Mix- 

ture 

ratio, 

(ft 

LOX flow rate, 
*0 

Com- 

bus- 

tion 

effi- 

ciency, 

V 

Fuel 

pressure - 
drop 
ratio, 
AP h 

P c 

LOX 

pressure - 
drop 
ratio, 

AP Io 

P C 

Lowe 

frequ 

clas 

cat 

r 

lency 

sifi- 

ion 

HigJ 

frequ 

clas 

cat 

tier 

ency 

sifi- 

Lon 

Ob- 

served 

lower 

fre- 

quency, 

Hz 

Ob- 

served 

higher 

fre- 

quency, 

Hz 

Fuel 

gain 

normal- 

izing 

factor, 

(* f ) 

1 ( ft - 5 

'“ o ’* 

LOX 

gain 

normal- 

izing 

factor, 

Normal- 
ized fuel 
pressure- 
drop 
ratio, 

r«\ 

Normal- 
ized LOX 
pressure- 
drop 
ratio, 

( AP *>\ 

lbj/in . 2 

N/m 2 

lb m /sec 

kg/sec 









v p *L 

v p c L 

74 

700 

4. 83 X 10 6 

5.01 

0. 561 

0.254 

0.787 

0.511 

0.656 

Stable 

Stable 

78 

195 

0.997 

1.00 

0.509 

0.656 

75 

696 

4. 80 

4.85 

.558 

.253 

.778 

.530 

.567 





78 

195 

1.06 

.973 

.562 

.551 

76 

696 

4.69 

4.85 

.567 

.257 

.748 

.577 

.478 





77 

195 

1.06 

.973 

.612 

.465 

78 

676 

4. 66 

4.77 

.573 

.260 

.732 

.287 

.942 





70 

205 

1.08 

.962 

.310 

.906 

80 

695 

4.75 

5. 16 

.583 

.264 

.760 

.248 

.576 





70 

190 

.950 

1.03 

.236 

.593 

82 

665 

4. 59 

4.76 

.538 

.244 

.766 

.184 

1.02 





80 

195 

1. 09 

.960 

.201 

.979 

83 

695 

4.75 

5.00 

.575 

.261 

.762 

.159 

.576 




80 

190 

1.00 

1. 00 

.159 

.576 

84 

675 

4.65 

4.80 

.567 

.257 

.740 

.182 

.622 




85 

195 

1.08 

.965 

.197 

.600 

85 

685 

4.72 

4.86 

.579 

.263 

.739 

.172 

.472 




80 

205 

1.05 

.975 

.181 

.460 

87 

695 

4.75 

4.70 

.555 

.252 

.773 

.133 

.834 




84 

165 

1.12 

.950 

.149 

.792 

88 

695 

4.75 

5. 17 

.585 

.265 

.758 

.122 

.522 


j 


83 

194 

.948 

1.03 

.116 

.538 

91 

693 

4.78 

4.97 

.591 

.268 

.737 

.0820 

.840 



i 

| 

82 

210 

1.01 

.995 

.0828 

.836 

92 

698 

4.81 

4. 80 

i .585 

.265 

.752 

.0850 

.440 

! 


! 

82 

195 

1.00 

1. 00 

.0850 

.440 

94 

680 

4.69 

5. 04 

.585 

.265 

.734 

.0480 

.831 


. 


-- 

--- 

.988 

1.01 

.0474 

.839 

96 

655 

4. 52 

’ 5.04 

.585 

.265 

.708 

.0760 

.386 




80 

200 

.988 

1. 01 

.0751 

.390 

98 

717 

4.94 

5.43 

.603 

.273 

.772 

| .0250 

1 .726 







.894 

1.07 

.0224 

.777 

100 

640 

4.41 

4.85 

.563 

.255 

.709 

I 

. 0250 

.430 



; 

79 

200 

, 1.06 

.973 

.0456 

.418 

108 

675 

4.65 

5.01 

.521 

.236 

I .817 

.711 

1.02 




: 63 

162 

.998 

1.00 

.0710 

1. 02 

109 

698 

4.81 

5.20 

j * 557 

.253 

, .801 

.698 

.532 



i 

78 

195 

.942 

1.04 

.0658 

.553 

111 

675 

4.65 

5.01 

.561 

.254 

> .759 

.756 

.444 


i 

! 

70 

190 

.998 

■ 1.00 

.0754 

.444 

113 

685 

4.72 

5.38 

.597 

.271 

■ .742 

.167 

.960 




78 

155 

.900 

1.06 

.150 

1.02 

114 

662 

4. 56 

4.71 

.542 

.246 

.755 

.196 

.960 




78 

200 

1.11 

.951 

.218 

.654 

115 

665 

4. 59 

5.13 

.585 

.265 

1 .723 

.196 

.688 




80 

198 

.959 

1. 03 

.188 

.480 

116 

655 

4.52 

4.88 

.561 

.254 

.730 

.214 

.397 


i 


75 

188 

1.05 

.980 

.225 

.389 

117 

617 

4.25 

4.81 

.563 

.255 

.682 

.240 

.361 



T 

71 

186 

1.07 

.967 

.257 

.349 



119 

687 

4.74x10*’ 

5. 01 

0. 567 

0.257 

0.764 

0.266 

0. 594 

Stable 

Stable 

72 

200 

0. 997 

1.00 

0.265 

0. 594 

123 

645 

4. 45 

5.30 

.567 

.257 

.732 

.081 

1.08 


68 

165 

.919 

1.05 

.0744 

1. 13 

124 

635 

4.38 

5.30 

.595 

.270 

.687 

.098 

.677 


78 

195 

.919 

1. 05 

.0900 

.711 

125 

592 

4. 08 

5.06 

.561 

.254 

.668 

.102 

.444 


80 

183 

.981 

1. 01 

.100 

.448 

126 

600 

4. 14 

4.83 

.536 

.243 

.697 

.112 

.362 

i 

80 

170 

1.06 

.970 

.119 

.351 

127 

565 

3.90 

4.73 

.530 

.240 

.660 

.115 

.327 


72 

175 

1. 10 

.955 

.127 

.312 

130 

705 

4.86 

5. 16 

.573 

.260 

.784 

.121 

.553 


80 

190 

.950 

1. 03 

.166 

.570 

131 

655 

4. 52 

5.01 

.561 

.254 

.736 

.166 

.431 


78 

190 

.997 

1.00 

.115 

.431 

135 

695 

4.75 

5.14 

.581 

.263 

.761 

.446 

.867 


78 

205 

.956 

1. 03 

.426 

.893 

136 

680 

4.69 

5.04 

.585 

.265 

.734 

.480 

.566 


74 

200 

.988 

1.01 

.474 

.572 

142 

618 

4.26 

4.81 

.558 

.253 

.689 

.080 

.384 


80 

195 

1.07 

.967 

.0856 

.371 

143 

596 

4.11 

4. 86 

.549 

.249 

.678 

.082 

.384 . 


82 

195 

1.05 

.975 

.0851 

.374 

144 

597 

4. 12 

5.01 

.542 

.246 

.695 

.080 

.284 


85 

182 

.997 

1.00 

.0798 

.284 

90 

635 

4.38 

5.10 

.597 

.271 

.675 

.422 

.207 

Unstable 

Unstable 

78 

150 

.969 

1.02 

.431 

.200 

97 

675 

4.65 

4. 87 

.585 

.265 

.721 

.053 

.218 

St at 

?le 

78 

198 

1.05 

.977 

.0556 

.213 

128 

515 

3.55 

4.88 

.542 

.246 

.594 

.136 

.209 


72 

178 

1.05 

.980 

.143 

.205 

145 

535 

3.69 

4.91 

.536 

.243 

.626 

.094 

.178 


72 

180 

.905 

.986 

.0861 

.176 

148 

487 

3.36 

4.96 

.536 

.243 

.571 

.078 

.093 


' 75 

178 

1.01 

.994 

.0788 

.0924 

71 

630 

4.34 

4.81 

.573 

.260 

.684 

.893 

.298 

Unstable 

Stable 

68 

202 

1.07 

.967 

.955 

.288 

77 

595 

4. 10 

5.15 

.618 

.280 

.613 

.744 

.324 


87 

205 

.953 

1.03 

.709 

.334 

81 

653 

4. 50 

5.03 

.610 

.277 

.676 

.321 

.406 


75 

210 

.990 

1.01 

.318 

.410 

89 

690 

4.76 

5. 19 

.610 

.254 

.785 

.152 

.435 


78 

153 

.943 

1. 04 

.143 

.452 

110 

623 

4.30 

4.92 

.585 

.265 

.667 

.895 

.396 


, 68 

200 

1.03 

.987 

.922 

.391 

112 

615 

4.24 

5. 12 

.558 

.253 

.700 

.927 

.390 



' 68 

200 

.962 

1.02 

.892 

.398 

121 

635 

4.38 

5.00 

.585 

.265 

.684 

.346 

.394 





70 

190 

1.00 

1. 00 

.346 

.394 

132 

598 



5.04 

.585 

.265 

.646 

.209 

.417 





74 

150 

.988 

1. 01 

.206 

.421 

133 

575 

— 

4.94 

.563 

.255 

.641 

.226 

.388 





71 

140 

1.02 

.990 

.231 

.884 

134 

535 

— 

5.01 

.561 

.254 

.601 

.214 

.276 





69 

205 

.998 

1. 00 

.214 

.276 

137 

650 

— 

5. 16 

.567 

.257 

.730 

.523 

.458 





74 

150 

.950 

1.03 

.497 

.472 

138 

615 

— 

4.77 

.540 

.245 

.707 

.573 

.439 





72 

147 

1.08 

.962 

.619 

.422 

139 

595 


5. 00 

.581 

.263 

.646 

.613 

.387 





70 

210 

1.00 

1.00 

.613 

.387 


to 

CJ1 



DO 

OS 


TABLE m. - CONFIGURATION B DATA 


Run 

Chamber pressure, 
P c 

Mix- 

ture 

ratio, 

<R 

i 

LOX flow rate, 

Com- 
bus- 
tion 
effl- 
ciency, 
: V c * 

Fuel 

pressure- 

drop 

ratio, 

AP h 

P c 

LOX 

pressure - 
drop 
ratio, 

* P Io 

P c 

Low 

frequ 

class 

cati 

er 

ency 

ifi- 

on 

Higher 

frequency 

classifi- 
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TABLE IV. - ESTIMATION OF COMBUSTION DELAYS FOR 


ANALOG COMPUTER SIMULATION 
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Figure 1. - Comparison of stability limit models for bipropellant rocket engines. 
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Figure 2. - Typical stability boundaries. Double dead-time model: fuel 
vaporization time delay, 0 millisecond; oxidizer vaporization time de- 
lay, 1.75 milliseconds; mixing time delay, 1,0 millisecond; theoreti- 
cal gas residence time, 0.7 millisecond; single dead-time model: time 
delay, 2.0 milliseconds; theoretical gas residence time, 0.7 millisec- 
onds. 
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Figure 3. - Rxket combustion laboratory flow system and test engine configuration. 
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Figure 4. - Cutaway diagram of hydrogen -oxygen chugging injector. 
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Figure 6. - Cold-flow fuel injector pressure-drop 
characteristics. 
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Figure 7. - Estimated fuel injector pressure-drop characteristics 
based on cold-flow tests. Fuel mass flow rate, 0. 11 pound mass 
per second. 
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Figure 9, - Typical observed chamber-pressure oscillation 
during chugging. Frequency, 42 hertz. 
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(a) Combustion-chamber static pressure, 650 psia (4.48xl0 6 
N/rrr); characteristic length, 91.0 inches (2.31 m) ; ob- 
served frequencies, 68 to 78 and 150 to 198 hertz. 



(b) Combustion-chamber static pressure, 300 psia (2.07xl0 6 
N/m 2 ); characteristic length, 95.5 inches (2.43 m); ob- 
served frequencies, 41 to 46 and 100 to 112 hertz. 


Figure 10. - Experimentally determined stability boundary. 
/Mixture ratio, 5.0; oxidizer mass flow rate, 0.55 pound 
mass per second (0.249 kg/sec). 
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Figure 11. - Determination of length required to vaporize 50 percent 
of oxidizer droplet mass. Combustion-chamber static pressure, 

300 psia (2.07x# N/m 2 ); mixture ratio, 5.0, ratio of cross- 
sectional area of combustion chamber to that of exhaust nozzle 
thrust, 16.8; combustion efficiency at 50 percent, 0.616; chamber 
length at 50 percent, 5.5 inches (14.0 cm); oxidizer jet velocity, 
735 inches per second (18. 7 cm/sec); time delay at 50 percent and 
oxidizer vaporization time delay, 7.48 milliseconds. 
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Figure 12. - Comparison of amplitudes of oscillation for experimental 
and simulation runs. Fuel injector pressure-drop ratio, 0.500, 
selected noise level. 
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Figure 13. - Effect of noise level on simulation 
stability boundary. Combustion-chamber 
static pressure, 650 psia (4.48xl0 6 N/m 2 ); 
fuel vaporization time delay, 0, oxidizer va- 
porization time delay, 4.4 milliseconds; mix- 
ing time delay, 2. 2 milliseconds; theoretical 
gas residence time, 2.0 milliseconds. 
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Fuel injector pressure-drop 



LOX injector pressure-drop ratio, APj Q /P c 


(a) Combustion-chamber static 
pressure, 650 psia (4.48xl0 6 
N/m 2 ); oxidizer vaporization 
time delay, 4.4 milliseconds*, 
mixing time delay, 2.2 milli- 
seconds; theoretical gas res- 
dencetime, 2.0 milliseconds. 


(b) Combustion-chamber static 
pressure, 300 psia (2. 07xl0 6 
N/m 2 ); oxidizer vaporization 
time delay, 8.0 milliseconds*, 
mixing time delay, 4.0 milli- 
seconds; theoretical gas resi- 
dencetime, 2. 1 milliseconds. 


Figure 14. - Comparison of experimental and simulation stability bound- 
aries. Fuel vaporization time delay, 0; noise level, 7 percent. 
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